Entanglement and purity of two— mode Gaussian states in noisy channels 
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We study the evolution of purity, entanglement and total correlations of general two-mode contin- 
uous variable Gaussian states in arbitrary uncorrelated Gaussian environments. The time evolution 
of purity, Von Neumann entropy, logarithmic negativity and mutual information is analyzed for a 
wide range of initial conditions. In general, we find that a local squeezing of the bath leads to a 
faster degradation of purity and entanglement, while it can help to preserve the mutual information 
between the modes. 
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I. INTRODUCTION 

In recent years, it has been increasingly realized that 
Gaussian states and Gaussian channels are essential in- 
gredients of continuous variable quantum information p| . 
Indeed, entangled Gaussian states have been successfully 
exploited in realizations of quantum key distribution 2] 
and teleportation protocols. 

In such experimental settings, the entanglement of a 
bipartite state is usually distilled locally, and then dis- 
tributed over space, letting the entangled subsystems 
evolve independently and move to separated spatial re- 
gions. In the course of this processes, interaction with the 
external environment is unavoidable and must be prop- 
erly understood. Therefore, the analysis of the evolu- 
tion of quantum correlations and decoherence of Gaus- 
sian states in noisy channels is of crucial interest, and has 
spurred several theoretical works 0, IE IE IE IE IE El • 

The evolution of fidelity of generic bosonic fields in 
noisy channels has been addressed in Ref. Q. In- 
deed, the relevant instance of initial two-mode squeezed 
vacua (possessing nontrivial entanglement properties) 
has drawn most of the attention in the field. Decoherence 
and entanglement degradation of such states in thermal 
baths have been analyzed in Refs. [E , whereas phase 
damping and the effects of squeezed reservoirs are dealt 
with in Refs. [EE US- In Ref. the author studies the 
evolution of a two-mode squeezed vacuum in a common 
bath endowed with cross correlations and asymptotic en- 
tanglement. Decoherence and entanglement degradation 
in continuous variable systems have been experimentally 
investigated in Ref. [Tl| . 

In this paper we address the general case of an arbi- 
trary two-mode Gaussian state dissipating in arbitrary 
local Gaussian environments. The resulting dynamics 
is governed by a two-mode master equation describing 
losses and thermal hopping in presence of (local) non 



classical fluctuations of the environment. 

We study the evolution of quantum and total corre- 
lations and the behavior of decoherence in noisy chan- 
nels. Quantum and total correlations of a state will be 
quantified by, respectively, its logarithmic negativity [l2T ] 
and its mutual information, while the rates of decoher- 
ence will be determined by following the evolution of the 
purity (conjugate to the linear entropy) and of the Von 
Neumann entropy. We present explicit analytic results, 
as well as numerical studies, on the optimization of the 
relevant physical quantities along the non-unitary evo- 
lution. Our analysis provides an answer to the ques- 
tion whether possible effective schemes to mimic general 
Gaussian environments 0, 0] are able to delay the de- 
cay of quantum coherence and correlations. We men- 
tion that, among such schemes, the most interesting for 
applications to bosonic fields is based on quantum non 
demolition (QND) measurements and feedback dynam- 
ics ^E 113 ■ We finally remark that the optimization 
of the quantities we are going to study with respect to 
phenomenological parameters turns out to be particu- 
larly relevant at 'small times', before decoherence has 
irreversibly corrupted the quantum features of the state, 
crucial for applications in quantum information. 

This paper is structured as follows. In Section II we 
provide a self-contained description of the general struc- 
ture of two-mode Gaussian states, including the charac- 
terization of their mixedness and entanglement. In Sec- 
tion III we review the evolution of Gaussian states in 
general Gaussian environments. In Section IV we focus 
on the evolution of purity and entanglement, determining 
the optimal regimes that can help preserving these quan- 
tities from environmental corruption. Finally, in Section 
V we summarize our results and discuss some outlook on 
future research. 



II. TWO-MODE GAUSSIAN STATES: 
GENERAL PROPERTIES 
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Let us consider a two-mode continuous variable sys- 
tem, described by an Hilbert space H = Hi ® H2 re- 
sulting from the tensor product of the Fock spaces Hi's. 
We denote by dj the annihilation operator acting on the 
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space Hi, and by Xi = (ai+a\)/y/2 &ndpi = (ai—a\)/V2 
the quadrature phase operators related to the mode i of 
the field. The corresponding phase space variables will 
be denoted by X{ and p,-. 

The set of Gaussian states is, by definition, the set of 
states with Gaussian characteristic functions and quasi- 
probability distributions. Therefore a Gaussian state is 
completely characterized by its first and second statisti- 
cal moments, which will be denoted, respectively, by the 
vector of first moments X = ((ah), (pi), (£2), (p2)) and 
by the covariance matrix er 



(£i) (Xj 



(1) 



First moments can be arbitrarily adjusted by local uni- 
tary operations, which do not affect any quantity related 
to entanglement or mixedness. Moreover, as we will see 
in Sec. IIHI they do not influence the evolution of second 
moments in the instances we will deal with. Therefore 
they will be unimportant to our aims and we will set 
them to in the following, without any loss of generality 
for our subsequent results. Throughout the paper, er will 
stand both for the covariance matrix and the Gaussian 
state g itself. 

It is convenient to express er in terms of the three 2x2 
matrices a, (3, 7 
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Positivity of g and the canonical commutation relations 
impose the following constraint for er to be a bona fide 
covariance matrix 16] 

a + l -U > , (3) 
where is the standard symplectic form 



n = 



Inequality Q is a useful and elegant way to express 
Heisenberg uncertainty principle. 

In the following, we will make use of the Wigner quasi- 
probability representation W, defined as the Fourier 
transform of the symmetrically ordered characteristic 
function [l7j . In Wigner phase space picture, the ten- 
sor product H = Hi ® H2 of the Hilbert spaces Hi's of 
the two modes results in the direct sum F = Ti ffi T 2 of 
the associated phase spaces iys. A symplectic transfor- 
mation acting on the global phase space T corresponds 
to a unitary operator acting on the global Hilbert space 
H • I n what follows we will refer to a transformation 
Si = Si © S2, with each Si S Sp(%m acting on Ti, as to 
a "local symplectic operation" . The corresponding uni- 
tary transformation is the "local unitary transformation" 
Ui = XJ\ ® U2, with each Ui acting on Hi- 

The Wigner function of a Gaussian state can be written 
as follows in terms of phase space quadrature variables 



W(X) = 



'X 1 



(4) 



where X stands for the vector (x\,pi,X2,P2) € T. 

It is well known that for any covariance matrix cr there 
exists a local canonical operation Si = Si © S2 which 
transforms <r to the so called standard form a s f [H^ ] 



STcrSi 
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States whose standard form fulfills a = b are said to be 
symmetric. Let us recall that any pure state is symmet- 
ric and fulfills c\ = — c 2 = \/ a? — 1/4. The correlations 
a, 6, Ci, and c 2 are determined by the four local sym- 
plectic invariants Deter = (ab — cf)(ab — c 2 ), Dcta: = a 2 , 
Det/3 = 6 2 , Det7 = cic 2 . Therefore, the standard form 
corresponding to any covariance matrix is unique. 

Inequality can be recast as a constraint on the 
Spurn invariants Deter and A(er) = Dcto; + Det/3 + 
2Det7: 



A(er) < - +4 Deter 



(6) 



Finally, let us recall that a centered two-mode Gaus- 
sian state can always be written as [2^, 



er = S T vS , 



(7) 



where S € Spum and v is the tensor product of thermal 
states with covariance matrix 



V = diag(n_, n_, n +1 n + ) . 



(8) 



The quantities n T form the symplectic spectrum of the 
covariance matrix er. They can be easily computed in 
terms of the Sp^^ invariants 



2ni 



A(er) =p \/A((t) 2 -4 Deter. 



(9) 



The symplectic eigenvalues encode essential informa- 
tions about the Gaussian state er and provide powerful, 
simple ways to express its fundamental properties. For 
instance, the Heisenberg uncertainty relation can be 
recast in the compact, equivalent form 

n->\. (10) 

A relevant subclass of Gaussian states we will make 
use of is constituted by the two-mode squeezed thermal 
states. Let S r — exp(|raia 2 — ^ra\a 2 ) be the two mode 
squeezing operator with real squeezing parameter r, and 
let — 1/(2 V / /I)1 be the tensor product of identical 
thermal states, where \i = Tr (p 2 ) is the purity of the 
state. Then, for a two-mode squeezed thermal state £ M r 
we can write ^ r = SrV^S^. The covariance matrix of 
£ M r is a symmetric standard form satisfying 



cosh 2r 
2VM 



ci 



-c 2 



sinh 2r 
2VM 



(11) 



and in the instance fi — 1 one recovers the pure two-mode 
squeezed vacuum states. Two-mode squeezed states are 
endowed with remarkable properties related to entangle- 
ment l22l l28j | ; their dynamics in noisy channels will be 
analyzed in detail. 
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Characterization of mixedness 



B. Characterization of entanglement 



Let us briefly recall that the degree of purity of a quan- 
tum state can be properly characterized either by the Von 
Neumann entropy Sv or by the linear entropy Si. Such 
quantities are defined as follows for continuous variable 
systems 



S v = - Tr(p lag) , 
Si = 1 - Tr(g 2 ) = 1 - /i 



(12) 
(13) 



where the purity \i = Tr(g 2 ) has already been intro- 
duced . We first point out that [i can be easily computed 
for Gaussian states. In fact, in the Wigner phase space 
picture the trace of a product of operators corresponds 
to the integral of the product of their Wigner represen- 
tations (when existing) over the whole phase space. Be- 
cause the representation of a state g is just W, for an 
n-mode Gaussian state we have, taking into account the 
proper normalization factor, 



TV 



W 2 d n xd n p 



1 



2™ V Deter 



(14) 



For Gaussian states, the Von Neumann entropy can be 
computed as well, determining their symplectic spectra. 
For single-mode Gaussian states, one has [24j 



Sv{<t) 



1 — n ( 1 + fi 



2^ 



In 



l-/x 



-In 



2[i 
1 + M 



(15) 



where \i can be computed from Eq. (|14|) for n — 1 . Sv is 
in this case an increasing function of the linear entropy, 
so that both quantities provide the same characteriza- 
tion of mixedness. This is no longer true for two-modes 
Gaussian states: in this case the Von Neumann entropy 
reads El EH 



S v (<t) = f[ii-(a)] + f[h+(a)] , 



(16) 



where 



f(x) = (x + I) \a{x + l)-( x -l) \n(x ~) 

and the symplectic eigenvalues n^er) are given by 
Eq. ©. 

Knowledge of the Von Neumann entropy Sv allows 
for the determination of the mutual information / de- 
fined, for a general bipartite quantum state g, as 1(g) = 
Sv(Qi) + Sv(@2) — Sv(q), where ft refers to the reduced 
state obtained tracing over the variables of subsystem 
j =/= i. The mutual information I(cr) of a two-mode 
Gaussian state er reads |2ll 



I(tr) = f(a)+f(b)-f(n-)-f{n+) 



(17) 



One can make use of such a quantity to estimate the 
amount of total (quantum plus classical) correlations con- 
tained in a state er |25| . 



We now review some properties of entanglement for 
two-mode Gaussian states. The necessary and sufficient 
separability criterion for such states is positivity of the 
partially transposed state er ("PPT criterion") [lfj. It 
can be easily seen from the definition of W(X) that the 
action of partial transposition amounts, in phase space, 
to a mirror reflection of one of the four canonical vari- 
ables. In terms of <Sp2,R© Sp2,R invariants, this results in 
flipping the sign of Det7. Therefore the invariant A(er) 
is changed into A(er) = A(cr) = Det a+ Det (3-2 Det 7. 
Now, the symplectic eigenvalues of er read 



\ 



A(er) =p A/A(er) 2 -4 Deter 



(18) 



The PPT criterion then reduces to a simple inequality 
that must be satisfied by the smallest symplectic eigen- 
value n_ of the partially transposed state 



1 

^2' 



which is equivalent to 



A(er) < 4 Deter + - 

4 



(19) 



(20) 



The above inequalities imply Det 7 = c\C2 < as a nec- 
essary condition for a two-mode Gaussian state to be 
entangled. The quantity n_ encodes all the qualitative 
characterization of the entanglement for arbitrary (pure 
or mixed) two-modes Gaussian states. Note that n_ 
takes a particularly simple form for entangled symmetric 
states, whose standard form has a = b 



n- = y/(a- \ci\)(a - |c 2 |) 



(21) 



As for the quantification of entanglement, no fully sat- 
isfactory measure is known at present for arbitrary mixed 
two-mode Gaussian states. However, a quantification of 
entanglement which can be computed for general two- 
mode Gaussian states is provided by the negativity Af, 
introduced by Vidal and Werner for continuous variable 
systems |l2j . The negativity of a quantum state g is de- 
fined as 



lelli-i 



(22) 



where g is the partially transposed density matrix and 
||d|| 1 = Tr V 6t<5 stands for the trace norm of an operator 
6. The quantity Af(g) is equal to | J^i -M> the modulus 
of the sum of the negative eigenvalues of g, and it quan- 
tifies the extent to which g fails to be positive. Strictly 
related to Af is the logarithmic negativity £jv", defined 
as Ej^f = ln||g||i. The negativity has been proved to be 
convex and monotone under LOCC (local operations and 
classical communications) [2(|, but fails to be continu- 
ous in trace norm on infinite dimensional Hilbcrt spaces. 
Anyway, this problem can be somehow eluded by restrict- 
ing to states with finite mean energy [28j • For two-mode 
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Gaussian states it can be easily shown that the nega- 
tivity is a simple function of n_, which is thus itself an 
(increasing) entanglement monotone; one has in fact [T^| 



En((t) = max{0, — ln2n_} 



(23) 



This is a decreasing function of the smallest partially 
transposed symplectic eigenvalue quantifying the 
amount by which Inequality l|19|) is violated. Thus, for 
our aims, the eigenvalue n_ completely qualifies and 
quantifies the quantum entanglement of a two-mode 
Gaussian state <x. 

We finally mention that, as far as symmetric states are 
concerned, another measure of entanglement, the entan- 
glement of formation Ep |29j. can be actually computed 
[23j. Fortunately, since Ep turns out to be, again, a de- 
creasing function of n_ , it provides for symmetric states a 
quantification of entanglement fully equivalent to the one 
provided by the logarithmic negativity EV- Therefore, 
from now on, we will adopt E_\f(cr) as the entanglement 
measure of Gaussian states, recalling that this quantity 
constitutes an upper bound to the distillable entangle- 
ment of quantum states [l2| . 



III. EVOLUTION IN GENERAL GAUSSIAN 
ENVIRONMENTS 

We now consider the local evolution of an arbitrary 
two-mode Gaussian state in noisy channels, in the pres- 
ence of arbitrarily squeezed ( "phase-sensitive" ) environ- 
ments. In general, the two channels related to the two 
different modes could be different from one another, each 
mode evolving independently in its channel. We will re- 
fer to the channel (bath) in which mode i evolves as to 
channel (bath) i. The system is governed, in interaction 
picture, by the following master equation [30| 

Q= E ~^iL[a]]0+^(iV' i + l)LMf ? 

i=l,2 

- | (MlD[a t ] g + Mi D[aJ]e) , (24) 

where the dot stands for time-derivative and the Lind- 
blad superoperators are defined by L[0}g = 20gO' [ — 
O^Og - gO^O and D[0]g = 20gO - OOg - gOO. The 
complex parameter Mi is the correlation function of bath 
i; it is usually referred to as the "squeezing" of the 
bath. Ni is instead a phenomenological parameter re- 
lated to the purity of the asymptotic stationary state. 
Positivity of the density matrix imposes the constraint 
\Mi\ 2 < Ni(Ni + 1). At thermal equilibrium, i.e. for 
Mi = 0, the parameter Ni coincides with the average 
number of thermal photons in bath i. 

A squeezed environment, leading to the master equa- 
tion Q24JI. may be modeled as the interaction with a 
bath of oscillators excited in squeezed thermal states 
[3l| . Several effective realizations o f sq ueezed baths have 
been proposed in recent years [l3l Il4|. In particular, in 
Ref. |l3T ] the authors show that a squeezed environment 
can be obtained, for a mode of the radiation field, by 



means of feedback schemes relying on QND 'intracavity' 
measurements, capable of affecting the master equation 
of the system. More specifically, an effective squeezed 
reservoir is shown to be the result of a continuous homo- 
dyne monitoring of a field quadrature, with the addition 
of a feedback driving term, coupling the homodyne out- 
put current with another field quadrature of the mode. 

Let gi = S(r.i, <pi)VfiiS(ri, <piY_be the environmen- 
tal Gaussian state of mode i |32|. Here Hi denotes 
the mean number of photons in the thermal state v^. 
Its knowledge allows to determine the purity of the 
state via the relation Hi = l/(2hi + 1). The operator 



jre l2v a t2 ) is the one-mode 



S(r, ip) = exp (ire 
squeezing operator. A more convenient parametrization 
of the channel, endowed with a direct phenomenological 
interpretation, can be achieved by expressing Ni and Mi 
in terms of the three real variables fii, Ti and ipi [33| 



cosh(2ri) 



y/(2N + l) 2 - 4|Mi| 



l + 4tf\Mi 



tan(2</?j) = -tan(ArgMi) 



(25) 



(26) 



(27) 



Note that the Gaussian state of the environment in bath i 
coincides with the asymptotic state of mode i, the global 
asymptotic state being an uncorrelated product of the 
states gi's, irrespective of the initial state. 

With standard techniques, it can be shown that the 
master equation (|24l) corresponds to a Fokker-Planck 
equation for the Wigner function of the system |3(j ■ In 
compact notation, one has 

W(X, t) = - [d x -X T + d x (Too d x ] W(X, t) , (28) 
with dx = (d Xl , d Pl , d X2 , d P2 ) and with a diffusion matrix 



= Vic 



' <T2oo — 



(T la 






(T2a 



(29) 



resulting from the tensor product of the asymptotic 
Gaussian states er ioo 's, given by 



= 2 



N t + Re Mi 
Im Mi 



Im Mi 
N - Re Mi 



(30) 



For an initial Gaussian state of the form Eq. the 
Fokker-Planck equation l|28|l corresponds to a set of de- 
coupled equations for the second moments and can be 
easily solved. Note that the drift term always damps 
to the first statistical moments, and it may thus be 
neglected for our aims. The evolution in the bath pre- 
serves the Gaussian form of the initial condition and is 
described by the following equation for the covariance 
matrix |i Igia. l34l 



(1 



,-rt\ 



<r(0)e 



-rt 



(31) 



This is a simple Gaussian completely positive map, and 
er(i) satisfies the uncertainty relation Eq. © if and only 
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if the latter is satisfied by both cr^ and er . The compli- 
ance of (Too with inequality Eq. J3J is equivalent to the 
conditions |M,| < N^Ni + 1). 

It is easy to see that Eq. (|5T|) describes the evolution 
of an initial Gaussian state <r in an arbitrary Gaussian 
environment cr^, which can in general be different from 
that defined by Eq. i|29|) . It would be interesting to find 
systems whose dynamics could be effectively described 
by the dissipation in a correlated Gaussian environment 
(recall that the instance we are analyzing involves a com- 
pletely uncorrelated environment). Some perspectives in 
this direction, that lie outside the scopes of the present 
paper, could come from feedback and conditional mea- 
surement schemes. 

The initial Gaussian state is described, in general, by 
a set of ten covariances. To simplify the problem and 
to better point out the relevant features of the non- 
unitary evolution, we will choose an initial state already 
brought in standard form: <t = cr s f. With this choice 
the parametrization of the initial state is completely de- 
termined by the four parameters a, b, c\ and C2, defined 
in Eq. JSJ. This choice is not restrictive as far as the 
dynamics of purity and entanglement are concerned. In 
fact, let us consider the most general initial Gaussian 
state <j evolving in the most general Gaussian uncorre- 
lated environment ®iO , i oa . The state <x can always be 
put in standard form by means of some local transfor- 
mation Si = (SiSi. Under the same transformation, the 
state of the environment ©i<x' ioo remains uncorrelated, 
with er'ioo = Sfcr ioo Si. All the properties of entangle- 
ment and mixedness for the evolving state are invariant 
under local operations. Therefore, we can state that the 
evolution of the mixedness and of the entanglement of 
any initial Gaussian state a in any uncorrelated Gaus- 
sian environment cr^ is equivalent to the evolution of the 
initial state in standard form Sf erSi in the uncorrelated 
Gaussian environment SjcTooSi. 

Finally, to further simplify the dynamics of the state 
without loss of generality, we can set cpx = ( correspond- 
ing to Im Mi = 0) as a "reference choice" for phase space 
rotations. 

Quite obviously, the standard form of the state is not 
preserved in an arbitrary channel, as can be seen from 
Eqs. (EOJ and (22). 

IV. EVOLUTION OF MIXEDNESS AND 
ENTANGLEMENT 

Let us now consider the evolution of mixedness 
and entanglement of a generic state in standard form 
(parametrized by a, b, c\ and C2) in a generic channel 
(parametrized by fix, rx, ^2, ti and t^)- Knowledge of 
the exact evolution of the state in the channel, given by 
Eq. I|31|) . allows to apply the results reviewed in SeclTTlto 
keep track of the quantities fi(t), S'y(t), I{t) and E_\f(t) 
during the non-unitary evolution in the channel. How- 
ever the explicit dependence of such quantities on the 
nine parameters characterizing the initial state and the 
environment is quite involved. We provide the explicit 
expressions m App. |U They give a systematic recipe 
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FIG. 1: Time evolution of logarithmic negativity of a non 
symmetric Gaussian state with a = 2, b = 1, ci = 1, C2 = — 1 
in several non correlated environments. The solid line refers 
to the case fii = y,2 = 1/2, n — T2 = 0; the dashed line refers 
to the case fii = 1/2, {12 = 1/6, r\ — T2 = 0; the dot-dashed 
line refers to the case fii = 1/6, /12 ~ 1/2, n = T2 = 0; the 
dotted line refers to the case fl\= [li = 1/2, n = T% = 1. In 
all cases the squeezing angle ip2 = 0. All the plotted quantities 
are dimensionless. 



to compute the evolution of mixedness, correlations and 
entanglement for any given Gaussian state in standard 
form (and, therefore, for any Gaussian state). 

We now investigate the duration and robustness of en- 
tanglement during the evolution of the field modes in 
the channels. Let us consider an initial entangled state 
<r e evolving in the bath. Making use of the separability 
criterion Eq. Q20[l. one finds that the state <r e becomes 
separable at a certain time t if 

u e" 4rt + v e" 3rt + w e~ 2Tt + y e~ r * + z = . (32) 

The coefficients u, v, w, y and z are functions of the 
nine parameters characterizing the initial state and the 
channel (see App. IA")i . Eq. (|32l) is an algebraic equation 
of fourth degree in the unknown k = e _rt . The solution 
k ent of such an equation closest to one, and satisfying 
k en t < 1 can be found for any given initial entangled 
state. Its knowledge promptly leads to the determination 
of the "entanglement time" t ent of the initial state in 
the channel, defined as the time interval after which the 
initial state becomes separable 

tent = — p In k ent ■ (33) 

The results of the numerical analysis of the evolution 
of entanglement and mixedness for several initial states 
are reported in Figs. ^ through |SJ In general, one can 
see that, trivially, a less mixed environment better pre- 
serves both purity and entanglement by prolonging the 
entanglement time. More remarkably, Fig. ^ shows that 
local squeezing of the two uncorrelated channels does not 
help preserving quantum correlations between the evolv- 
ing modes. Moreover, as can be seen from Fig. ^ states 
with greater uncertainties on, say, mode 1 (a > b) better 
preserves its entanglement if bath 1 is more mixed than 
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FIG. 2: Time evolution of logarithmic negativity of a symmet- 
ric Gaussian state with a = 1.5, b — 1.5, c\ = 1.2, C2 = —1.4 
in several non correlated environments. The solid line refers 
to the case [i\ — fi2 = 1/2, n = r2 = 0; the dashed line refers 
to the case fj,i = fi2 = 1/4, n = r2 — 0; the dot-dashed line 
refers to the case fi\ = fi2 — 1/2, r\ — T2 — 1; the dotted 
Une refers to the case /Mi = [i2 = 1/2, n = 0, T2 = 1.5. In all 
cases the squeezing angle if 2 = 0. All the plotted quantities 
are dimensionless. 
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FIG. 4: Time evolution of logarithmic negativity of a two- 
mode squeezed thermal state with initial purity /1 = 1/9, 
r = 1 in several non correlated environments. The solid line 
refers to the case /Lti = fi2 = 1/2, n = r-2 = 0, </72 = 0; the 
dashed line refers to the case /ii = 4, ^2 = 1, ri = r2 — 0, 
tf2 = 0; the dot-dashed line refers to the case /ii = /i2 = 1/2, 
ri = r2 = 1, (^2 = 0; the dotted line refers to the case /ii = 
fi2 = 1/2, n = ?"2 = 1, 1^2 = 7r/4. All the plotted quantities 
are dimensionless. 
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FIG. 3: Time evolution of logarithmic negativity of a two- 
mode squeezed state with r = 1 in several non correlated 
environments. The solid line refers to the case fit = \i% = 1/2, 
ri = T2 = 0, if2 = 0; the dashed line refers to the case /ii = 4, 
M2 = 1, fi = T2 = 0, f2 = 0; the dot-dashed line refers to 
the case /ii = [12 = 1/2, n = r2 — 1, ^2 = 0; the dotted line 
refers to the case fii — /12 = 1/2, n = r% = 1, (^2 = 7r/4. All 
the plotted quantities are dimensionless. 

bath 2 (/ii < H2). A quite interesting feature is shown 
in Fig. |SJ the mutual information is better preserved in 
squeezed channels, especially at long times. This prop- 
erty has been tested as well on non entangled states, en- 
dowed only with classical correlations, see Fig. [HJ and on 
two-mode squeezed states, see Fig. [7] and seems to hold 
quite generally. In Fig. [3 the behavior of some initially 
symmetric states is considered. In this instance we can 
see that, in squeezed baths, the entanglement of the ini- 
tial state is better preserved if the squeezing of the two 
channels is balanced. 

The analytic optimization of the relevant quantities 
characterizing mixedness and correlations in the channel 
turn out to be difficult in the general case. Thus it is 
convenient to proceed with our analysis by dealing with 
particular instances of major phenomenological interest. 



FIG. 5: Time evolution of the purity of two-mode squeezed 
thermal states. The solid line refers to a two-mode squeezed 
vacuum state with r = 1 in an environment with /Lti = [12 = 
1/2, n = T2 — 0; the dashed line shows the behavior of 
the same state for fi\ — [12 — 1/2, ri = f2 = 1. The dot- 
dashed line refers to a mixed state with /i = 1/9, r — 1 in an 
environment with fii = fj,2 = 1/2, n = T2 — 0; the dotted line 
refers to the same state for /ii = fi2 — 1/2, n = f2 = 1. The 
squeezing angle if 2 has always been set to 0. All the plotted 
quantities are dimensionless. 



A. Standard form states in thermal channels 

In this subsection, we deal with the case of states in 
generic standard form (parametrized by a, 6, Cj) evolving 
in two thermal channels (parametrized by two - generally 
different - mean photon numbers Ni's). This instance is 
particularly relevant, because it gives a basic description 
of actual experimental settings involving, for instance, 
fiber-mediated communications. 

The purity fi of the global quantum state turns out to 
be a decreasing function of the Ni's at any given time. 
The symplectic eigenvalue n_ is also in general an in- 
creasing function of the Ni's. Therefore, the entangle- 
ment of the evolving state is optimal for ideal vacuum 
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FIG. 6: Time evolution of the Von Neumann entropy of sev- 
eral Gaussian states in a thermal environment with /ii = 
fj,2 = 1/3. the solid line refers to a state with a = 1, 
b = ci = — C2 = 1; the dashed line refers to a two-mode 
squeezed vacuum state with r = 1; the dot-dashed line refers 
to a squeezed thermal state with /j, = 1/16, r = 1; the dot- 
ted line refers to a non entangled state with a — b = 2, 
ci = — C2 = 1.5. The squeezing angle ip2 has always been 
set to 0. All the plotted quantities are dimensionless. 
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FIG. 7: Time evolution of the mutual information of two- 
mode squeezed thermal states in an environment with fii — 
|U2 = 1/3- The solid line refers to a pure state with r = 1 in a 
non squeezed environment; the dotted line refers to the same 
state in an environment with n = r% = 1; the dashed line 
refers to a squeezed thermal state with /j, — 1/16, r — 1 in a 
non squeezed environment; the dot-dashed line refers to the 
same state in a squeezed environment with n = ri = 1. The 
squeezing angle if2 has always been set to 0. All the plotted 
quantities are dimensionless. 



environments, which is quite trivial, recalling the well 
understood synergy between entanglement and purity for 
general quantum states. 



FIG. 8: Time evolution of the mutual information of Gaussian 
states in an environment with /Ui = fi2 = 1/3. The solid line 
refers to state with a = 2, b = ci = — C2 = 1 in a non squeezed 
environment; the dotted line refers to the same state in an 
environment with n = rz = 1; the dashed line refers to a 
non entangled state with a = b = 2, ci = — C2 = 1.5 in a 
non squeezed environment; the dot-dashed line refers to the 
same state in a squeezed environment with n = Ti = 1. The 
squeezing angle (p% has always been set to 0. All the plotted 
quantities are dimensionless. 



symmetric entangled state maintains its symmetric stan- 
dard form if evolving in equal, independent environments 
(with Ni = N2 = Nb)- This is the instance we will con- 
sider in the following. 

Let us suppose that |ci| < |c 2 |, then Eqs. (jTT?|l and ij?T|) 
provide the following bounds for the entanglement time 



In 1 



\ci\ - a 



1 



C2 



< Ttent < In i 

i\b 

(34) 

Imposing the additional property c\ — —c^ we obtain 
standard forms which can be written as squeezed thermal 
states (see Eqs.^^. For such states, Inequality Eq. 
reduces to 



tent — p hi 



1 



(35) 



In particular, for fj, = 1, one recovers the entanglement 
time of a two-mode squeezed vacuum state in a thermal 
channel 0, El E| • Note that two-mode squeezed vac- 
uum states encompass all the possible standard forms of 
pure Gaussian states. 



C. Two— mode squeezed thermal states 



B. Entanglement time of symmetric states 

We have already provided a way of computing the en- 
tanglement time of an arbitrary two-mode state in ar- 
bitrary channels. The expression of such a quantity is, 
unfortunately, rather involved in the general case. How- 
ever, focusing on symmetric states (which satisfy a = b), 
some simple analytic results can be found, thanks to the 
simple form taken by n_ for these states. An initially 



As we have seen, two-mode squeezed thermal 
states constitute a relevant class of Gaussian states, 
parametrized by their purity /1 and by the squeezing pa- 
rameter r according to Eqs. i|llfl . In particular, two- 
mode squeezed vacuum states (or twin-beams), which 
can be defined as squeezed thermal states with fj, = 1, 
correspond to maximally entangled symmetric states for 
fixed marginal purity. Therefore, they constitute a cru- 
cial resource for possible applications of Gaussian states 
in quantum information engineering. 
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For squeezed thermal states (chosen as initial condi- 
tions in the channel), it can be shown analytically that 
the partially transposed symplectic eigenvalue n_ is at 
any time an increasing function of the bath squeezing an- 
gle if2- "parallel" squeezing in the two channels optimizes 
the preservation of entanglement. Both in the instance of 
two equal squeezed baths (i.e. with r± = T2 = r) and of 
a thermal bath joined to a squeezed one (i.e. r\ = r and 
f2 =0), it can be shown that n_ is an increasing function 
of r. These analytical results agree with those provided 
in Ref. @ in the study of the qualitative degradation of 
entanglement for pure squeezed states. The proofs of the 
above statements are sketched in App. [BJ 

Such analytical considerations, supported by direct nu- 
merical analysis, clearly show that a local squeezing of 
the environment faster degrades the entanglement of the 
initial state. The same behavior occurs for purity. The 
time evolution of the logarithmic negativity of two-mode 
squeezed states - thermal and pure - is shown in Figs.|3] 
and 2] The evolution of the global purity is reported in 
Fig. [SI while the evolutions of the Von Neumann entropy 
and of the mutual information are shown, respectively, 
in Figs. and H 



in environments with lower average number of photons. 

The present study may be be extended to the case of 
n-mode Gaussian states. This generalization would be 
desirable, since the practical implementation of quantum 
information protocols usually requires some redundancy. 
For three-mode Gaussian states, separability conditions 
analogous to Inequality Eq. 10 have been determined 
|35j . and could be exploited to provide a qualitative pic- 
ture of the evolution of three-mode entanglement in noisy 
channels. 
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APPENDIX A: EXPLICIT DETERMINATION OF 
MIXEDNESS AND ENTANGLEMENT IN THE 
GENERAL CASE 



V. SUMMARY AND CONCLUSIONS 

We studied the evolution of mixedness, entanglement 
and mutual information of initial two-mode Gaussian 
states evolving in uncorrelated Gaussian environments. 
We derived exact general relations that allow to deter- 
mine the time evolution of such quantities, and provided 
analytical estimates on the entanglement time. The op- 
timal bath parameters for the preservation of quantum 
correlations and purity have been determined for thermal 
baths and for two-mode squeezed states in more general 
baths. A detailed numerical analysis has been performed 
for the most general cases. 

We found that, in general, a local squeezing of the 
baths does not help to preserve purity and quantum 
correlations of the evolving state, both at small times 
(i.e. for Ft < 1) and asymptotically. On the other hand, 
local squeezing of the baths can improve the preservation 
of the mutual information in uncorrelated channels. Be- 
sides, coherence and correlations are better maintained 



Here we provide explicit expressions which allow to 
determine the exact evolution in uncorrelated channels 
of a generic initial state in standard form. The relevant 
quantities £V, //, Sv, I, as we have seen, are all functions 
of the four Sp(2,R) ® Sp(2,R) invariants. Let us then write 
such quantities as follows 

4 

Deter = ^2^ k e~ krt , (Al) 

fc=0 
2 

Deta = ^a fe e- fcrt , (A2) 

k=0 
2 

Det/3 = 5> fc e- fcrt , (A3) 

fc=0 

Det 7 = 72 e" 2rt , (A4) 

(A5) 

defining the sets of coefficients £j, a t , ft, ji. One has 



2.2 



a z b 



+(c? 



b 2 2 cosh2r2 2 cosh2ri cosh 2ri cosh 2r2 cosh2ri cosh2r2 
a b ab h ab a — : ^ b- 



/'2 



Ml 



cosh2r2 6 cosh 2?^ cosh 2r*i cosh 2r 2 sinh 2r^ sinh 2r 2 cos 2<p 2 



2/l 2 



2/12 



2/*i 
+ b- 



4^iA*2 



4/il/X2 



— ab 



2 1 sinh 2r2 cos 2(^2 sinh2r*i sinh 2r\ cosh 2r2 cosh 2r\ sinh 2r2 cos 2(^2 
+ (ci - c 2 " ■ 



2//i 



4/i 1M2 



4^lM2 



■ 2 2 
+c x c 2 



1 



1 b 2 2 cosh2r2 2 cosh2ri 

' ab 



4M| 



4M? 



a z b- 



fJ-2 



Ml 



cosh 2ri cosh 2r2 cosh2ri cosh2r2 
2ab h 3a— ^ — H 3o- 



M1M2 



4/iiM2 



4/iiM2 



(A6) 
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Si 



-(c? + c§) 

-(c? + ci) 

-(c?-ci) 
cosh 2ri 

i 2~~ 

1 

16/z?m| ' 



sinh 2r2 cos 2^2 sinh 2r\ sinh 2ri cosh 2r2 
a 1- o 2- 



2^2 



2^1 



4^1 



) cosh 2ri sinh 2r2 cos 2^2 

4/Lti/i 2 



cosh2r2 6cosh2ri cosh 2ri cosh 2r2 sinh 2ri sinh 2r2 cos 2^2 



ab 



2fi 2 2^! 
cosh2ri cosh2r2 



4^iM2 



4/^1/^2 



3a 



cosh2ri 



36 



cosh 2r2 



M1M2 4/ii/i^ 4^J^ 2 

cosh 2ri cosh 2r2 sinh 2ri sinh 2r2 cos 2(f2 

4mi/x 2 4^2 

sinh 2ri cosh 2r 2 cosh 2r\ sinh 2r2 cos 2(f2 



+ b 



4/Ltl/Z 2 

cosh 2r 2 

4^?^2 4/^/x 2 



4^i ^2 



1 



I6M1M2 



2, ,2 ' 



4Mi M2 



(A7) 



(A8) 
(A9) 
(A10) 



Ot2 
ft 

72 



cosh 2r\ 1 

a !~ 1 — 2 

Mi 4/^ 



cosh2ri 



/'1 



1 



4/i? 



&2 & cosh2r 24 

M2 

, cosh2r2 1 
& - 2 



4m! 



/'2 



4^1 ' 
cic 2 . 



4^! 



(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 
(A17) 



u = 


s 4 


; 




V = 


S3 


; 




w = 


£2 


— «2 


-(3 2 - 


y = 


Si 


- ai 


- Pi , 


z = 


S 


- "0 


-0o + 



The coefficients of Eq. (|32|l , whose solution k en t allows 
to determine the entanglement time of an arbitrary two- 
mode Gaussian state, read 



(A18) 
(A19) 

I72I , (A20) 
(A21) 

\ ■ (A22) 



APPENDIX B: PROOFS FOR TWO-MODE 
SQUEEZED STATES 

In this appendix we consider a two-mode squeezed 
thermal state of the form of Eq. i|ll[) as the initial in- 
put in the noisy channels. 

We first deal with the dependence of entanglement and 
mixedness on the squeezing angle tf2 of bath 2. It can be 
easily shown (see App.[XJ that A(er) does not depend on 
(fi2, whereas Deter turns out to be a decreasing function 
of cos if2- Therefore, since the symplectic eigenvalue n_ 
increases with Det er, one has that if2 = is the optimal 



choice for maximizing both entanglement and purity of 
the evolving state. 

We now address the instance of two equally squeezed 
baths, with Ni = Nb, r% = tb and Lp 2 — 0. The time 
dependent covariance matrix er 2m can be written in the 
form 



& 2r. 




with 



k 



cosh 2r 
sinh 2r 



-11 



-ri 



+ (N B + -)e^(l- c - 



■rt\ 



The standard form of cr 2m is easily found just by squeez- 
ing the field in the two modes of the same quantity 
V 3+ 1 3- ■ The result is a symmetric standard form, whose 
smallest partially transposed symplectic eigenvalue n_ 
can be computed according to Eq. (EIJ 

ri_ = — k)(j + — k) = dcosh2rs + . . . , 

where the terms that do not depend on tb are irrelevant 
to our discussion and have thus been neglected. The 
coefficient d is a positive function oft, r and Nb, so that 
the best choice to maximize entanglement at any given 
time is given by rs = 0. Quite obviously, n_ turns out 
to be an increasing function of Nb as well. 

Finally, we deal with the instance in which bath 1 is 
squeezed while bath 2 is thermal, with r2 — 0. For ease 
of notation we define lerl = Deter. We recall that 2n 2 _ = 



A — y A 2 — 4|er|. Thus, for entangled states (for which 
n_ < 1/2), one finds 

d W \(2h 2 _) > -Ad K {2h 2 _) > . 
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The sign of the quantity 4<9 ri |er| — d ri A for the case of 
the initial two-mode squeezed can be shown, after some 
algebra, to be determined by 

4( e r * - 1) cosh 2rn\ + (3 + cosh 4r)n 2 - ( e r * + 1) cosh 2r . 



This second degree polynomial is positive for ri2 = N2 + 
1/2 > 1/2. This proves that the entanglement decreases 
as the squeezing of bath 1 increases. 
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